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Abstract 

The free energy of the supersymmetric t — J model is expressed in terms of finite tem- 
perature excitations above thermodynamic equilibrium. This reveals that the free energy 
has the form of noninteracting fermions with temperature dependant excitation spectra. 
We also discuss the ground state and zero temperature excitations. 

PACS numbers:75.10-b 71.20AD 



1 Introduction 



Strong electron correlations are believed to be important in understanding the high T c supercon- 
ductors [Q] . This idea is supported by the fact that the high T c compounds display antiferromag- 
netism in the absence of doping. The t — J model describes strongly correlated electrons with 
antiferromagnetic exchange interactions, and has been proposed as a candidate to model high 
T c superconductors by Zhang and Rice 0. The one dimensional version of the model becomes 
integrable || and supersymmetric [|j] for special values of the parameters t and J. Integrability 
allows the model to be solved by Bethe's ansatz and supersymmetry leads to the construction of 
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three distinct Bethe ansatz solutions. Lai || and Sutherland || each found a solution in the con- 
text of models of hard core bosons and fermions which are now understood to be equivalent to the 
super symmetric t — J model. Schlottman |J applied the Bethe ansatz to the t — J hamiltonian, 
and found Lai's solution describes the model when electrons are treated as dynamical objects in 
a background of empty sites. Sarkar |7| discovered that Sutherland's solution applies to the t — J 
model if one treats holes and spin down electrons as dynamical objects in a background of spin 
up electrons. There is a third and quite recent solution due to Essler and Korepin H, which is 
similar to Sutherland's, but interchanges the roles of the spin down electrons and the holes. 

The ground state of the t—J model was constructed in [||, ||, [7|]. The free energy was calculated 
in [[J. Bares, Blatter, and Ogata |J use both Lai's and Sutherland's ansatzes to give a detailed 
account of the ground state and excitation spectrum, including extensive numerical analysis. 

The outline of the paper is as follows: in sec. 2 we define the model and take the thermo- 
dynamic limit of the the Bethe ansatz equations. In sec. 3, the bulk free energy is expressed 
in terms of finite temperature excitations. In these terms, the free energy has the form of a 
noninteracting system with temperature dependant excitation spectra. This structure has been 
observed in such models as the 5-function bose gas JlO| and the XXZ spin chain |ll|], here these 
ideas are extended to th e t — J model. In sec. 4 we discuss the ground state and the order one 
excitation spectrum. The results agree with those of @, |J, we present them here for complete- 
ness, and to see how see how they arise from a different form of Bethe ansatz. In the appendix 
we note that the different ansatzes are related by an interchange of particle and hole rapidities. 



2 Formulation 

The t — J model describes spin-| fermions on a lattice. The on-site coulomb repulsion is taken 
to be infinite, so that there is no double occupancy of the sites. The constraint of no double 
occupancy implies that the number of empty sites (holes) is conserved, which allows one to think 
of the empty sites as dynamical objects. The parameter t describes how easily the electrons hop 
from site to site and J describes the strength of an antiferromagnetic exchange interaction. The 
Hamiltonian is 

H = V iJ2 £ + <Wfc,) + J(4 ■ S J+1 - \n n ]+1 )}V + 2N-L. (2.1) 

J = l <7=U 

L is the number of lattice sites. c J)(J , and Cy a , are electron creation and annihilation operators 
which obey the anti-commutation relations 

{Cj,a, Cj',o>} = 0, {C j)ff , ct, CT ,} = SjjiS^, <T, (J* =|, | . (2.2) 
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V is the operator which projects out doubly occupied sites, Sj = cij a a ai(7 >Cj ja ', where a a y is the 
vector of Pauli spin matrices, and rij is the number of electrons at site j. N is the operator for the 
total number of electrons. We have added the term 2N — L onto the definition of the Hamiltonian 
in H [J, which just shifts the energy and chemical potential. We shall consider this model at 
the special point J = It = 2, where the model is integrable and supersymmetric Also, at 
this special point the Hamiltonian (|2.1| ) may be rewritten as a graded permutation operator (?]], 

tf=-£lW, (2.3) 

i=i 

where interchanges the states on neighboring sites, with a minus sign if electrons sit on 

both sites. 

The eigenstates and spectrum of the Hamiltonian (|2.1|) at the supersymmetric point may 
be found by nested Bethe ansatz, a technique introduced by Yang in |T2[. The spectrum of 



the Hamiltonian (|2.1| ) and other conserved quantities are given in terms of a set of coupled 
algebraic equations known as Bethe ansatz equations(BAE). As mentioned in the introduction, 
there are three distinct sets of BAE for the model. We shall analyze the solution due to Essler 
and Korepin ||, where the coordinates of holes and spin down electrons are taken as dynamical 
objects moving in a background of spin up electrons. This solution expresses the eigenvalues of 
the Hamiltonian and total momentum as: 

Ni+N h , 

E = L- £ ^— (2.4) 
f=t v A/-* y 

where the Nh + spectral parameters A; describe the motion of holes and spin down electrons. 
They must satisfy the Bethe ansatz equations 

(^)'=fi^ 1 = 1 J^+AW, (2.6) 

1= n Y^ 11 >' ' V. (2.7) 
i=i k - A/3 - i 

Ni, iVt, and Nh — L — N± — are the number of spin up electrons, spin down electrons, and 
holes respectively. The iVj_ parameters describe the motion of the spin down electrons relative 
to the holes. 



3 



The equations ( |2.6| ) and (57?) have real and complex solutions. The complex solutions are 
of a special form known as strings, which may be found by fixing Nh and iVj_ and letting the 
lattice size L go to infinity. One finds complex solutions consisting of n A's and n — 1 A's in the 
following combinations ||: 



Xj = X + i(n + 1 — 2j) j — l,...,n 
A s = \ + i(n-25) 5 = 1, - 1, 



(2.8) 



for arbitrary n = 1, oo. The real part common to all the parameters, A, is known as the center 
of the string. Physically, string solutions are bound states in the sense that an eigenfunction 
with the string combinations decays exponentially with respect to the coordinates of the down 
spins and holes. We shall see later in this section that the string solutions (|2.8|) yield a complete 



spectrum of the hamiltonian. In a finite box, the string solutions ( |2.8|) are not exact. We assume 



along with many authors (ie. ||T3| , [14], [iTJ]), that the corrections vanish in the thermodynamic limit. 

We wish to write BAE for the centers of the strings. Let X™ the center of the ath string of 
length n, and be real rapidities for spin down electrons not in any bound state. Denote the 
number of strings of length n by N n and the number of free spin downs by N\. Multiplying 
together the equations for each part of the string yields 



•A" 



+ in\ l 



■X'X-in 



oo N n 

II II F nm(K 
m=l 7=1 

(772,-7)7^(71,0:) 



^^XZ-Ap-in 



(2.9) 



where n = 1, 00 and 



1 = tt n p - a 

An - X n - in ' 



where 



X 



X 



X 



\n — m\ 



\n — ml + 2' 



x • • • x e 



'n + m — T 



(2-10) 
(2.11) 



+1 



and we have defined e(x) = - 

The energy and momentum in terms of the centers of the strings are 

An 



OO N n 



and 



^ ^ (X n ) 2 + n 2 

n=l a=l \^a) ~ 11 



oo N n \n 



n=l a=l 



n 



(2-12) 



(2-13) 
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where p(X) = 2 tan l (X) — tt. 

In order to count the solutions and take the thermodynamic limit of the BAE, we take the 
logarithm of Q and (PHQ : 

\n oo Nm N A \ n A 

L9(±) = +EE <SW(A£ - A-) + E £), (2.14) 



n m =l 7 =i /3=i n 



and 



EE^)=H, (2-15) 



n=l " ra 



where #(:r) = 2 tan 1 (x), and 



f^(HS) + 2g( , nzn ) + • • • + 2fl(-r 3L - 5 ) ifn^m 
nml j ~ 1 26(1) + 20(f) + ■ • • + 26(^- 2 ) if n = m. ^ iDJ 

7™ and Jg are integers (half integers) arising from the choice of the branch of the logarithm. J™ 
is an integer (half integer) if (L — Y^m^n N m — N A ) is even (odd). Jp is an integer (half integer) 
if (En=i^n) is even (odd). 

Given a solution {A™}, {Ap} define functions 



*»(A) = d{-) - - E E ©-(A - A™) - y E *( ") (2-17) 

" ^ m=l 7=1 ^ /3=1 72 

for n = 1, oo, and 

i oo N„ a \n 

^ A (A) = }EE^ (2-18) 

^ „=1 a =i n 

By definition, Lz n (\ n a ) = 2nl% and Lz A {A l3 ) = 2irJ p . If z n {\) and 2 A (A) are monotonic functions, 
then specifying sets of integers and {Jp} uniquely determines the solution of Q2.9Q and 

( |2.10D . Define the integers corresponding to infinite spectral parameter by 2^1^ = Lz 11 (oo) 
and 27rJoo = L,2 A (oo). In order to count the states correctly, we discard solutions with infinite 
spectral parameter, yielding the maximum allowable integers 

jn = jn _ 1 

max oo 

J oo 

= -(L- E N m (t nm -l)-N A -2), (2.19) 

Z m=l 
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and 



i 



(2.20) 



1 oo 

5(E^»- 2 ). 

Z n=l 



where t nm = 2min(m, n) — 5 nm . Counting solutions of the BAE ( |2.9| ) and ( 2.10| ) using the 



above, and taking into account that each Bethe state generates a super symmetry multiplet of 
size d = 8S Z , one obtains a sum for the total number of states identical to the corresponding 
formula found by Foerster and Karowski using Sutherland's ansatz in [IB]. There it is shown 
that the sum adds up to 3 L , the total number of possible states. Thus there are no additional 
complex solutions other than (|2.8|) contributing to the spectrum. 



z n and z 



We may also use equations ( [2.171) an d (|2.18| ) to define hole rapidities. Having assumed that 
A are monotonic functions of A and A, then for any unoccupied integer I™, there exist a 
corresponding hole rapidity A™. We will call the combined sets of {A™} and {A™} the vacancies 
for parameters. Vacancy integers run through the entire allowed range of integers. 



2.1 Thermodynamic Limit of BAE 

We now take the thermodynamic limit of the Bethe ansatz equations (|2.9| ) and ( 2.10| ), fixing the 
electron and magnetization densities. The solutions to the BAE become densely packed with 
differences between neighboring A™ +1 — A™, Aj + i — Aj ~ 0(1/ L). One passes to a description of 
the states in terms of densities. Define particle and hole densities: 

p£(A) = lim 1 p h n {\) = lim 1 (2.21) 

L^OO i^Api +1 — A/'i ) L^OO Ll\AJn^ — Ajn ) 

and 

a?(A) = r 1 ™ TTK 1 A V Cj/l(A) = r lim 771 1 A V (2 ' 22) 

l^oo L(Aj. +1 - Aj.) l^oo L(A J]+1 - A j.) 

and the total density of vacancies is p l n = p v n + and a t = a p + a h . In the thermodynamic limit, 
the Bethe ansatz equations become integral relations between the densities: 

°° rao poo 

Pn(A) = /»(A) - E / d\'A n M ~ A')/4(A') - / dAf n (X - A)o*(A) (2.23) 



m=l 

and 



°° POO 

a t (A) = Y, d\f n (A-\)(? n (\), (2.24) 

n=l J -°° 
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where 

/»(*) = ^Jx 6 ^ = -WT-2 ( 2 - 25 ) 

Z7T d A 7T \ z + 7r 



and 



j / \ \ _ I f\n-m\\X) + 2/| n _ m | +2 + • • • + 2/„ +m _ 2 if 7i 7^ 771 ^ 2g") 

2/2 + 2/4 H h 2/ 2n _ 2 if n = m. 



3 Thermodynamics 

3.1 Free Energy 



To calculate the free energy, we follow the techniques developed by Yang and Yang [HJ, Taka- 
hashi ||13||, and Gaudin jTJJ. The free energy is a functional of the densities: 



F — E — fiN — BS Z — TS, (3.1) 

where the energy is (all integrals beyond this point are between — oo and oo unless explicitly 
otherwise) 

00 „ 

E = L(1-J2 rfA4vr/ n (AK(A)), (3.2) 

n=l 

the number of particles is 

~ OO r. 

N = L(1- dXY, P&(A) + / dAa?(A)), (3.3) 

n=l 

the magnetization is 

OO p „ 

2S Z = L(l - £(2n - 1) J d\f? n {\) - j dAa p (A)), (3.4) 



n=l 

and the entropy is given by 



S = L 



dX £ [pi hgipi) - (f n tog(pP) - p h n hg(p h n )} + / dA[a' log(a') - a" log(^) - a h \og(a h )] 



n=l 



(3.5) 

To obtain the entropy ( |3.5| ), note that the number of states in in an interval dX in one of the 
fermi seas is 
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This the number of ways to put [Lp p d\] "particles" into [Lp t dX\ vacancies. Using Stirling's 
approximations for the factorial function at large N,\ogN\ ~ NlogN, and adding up the con- 
tribution due to all types of particles, one arrives at (BT5[) . 



We now have the free energy for an arbitrary state. To find the equilibrium value, one 
minimizes the free energy with respect to the densities. We take the pP n and a p as the independent 
variables, with p l n and a 1 determined from ( |2.23| ) and ( |2.24| ). Setting the variation of the free 



energy equal to zero yields a infinite set of coupled nonlinear integral equations for the equilibrium 
densities: 

e n (X) = e° (A) + T £ j dvA nm {\ - u) log(l + e"^) —T J dAf n (\ - A) log(l + e~^), (3.7) 



m=l 

with n = 1 oo, and 



oo „ 

e A (A) = -p-B + Tj2 d\f n (A - A) log(l + e"^), (3.8) 

n=l 

where we have written the equations in terms of the functions 

e n = Tlog(4), eA^TlogA, (3.9) 
Pn a p 

and 

e° n (A) = -4vr/ n (A) - B{2n - 1) + p. (3.10) 
Evaluating the free energy at the minimum yields the equilibrium value 

F „ „^ r „ „ t „w , 

~L 



C (X) 

[l + B-p)-Tj2l dXfn(X) log(l + e"V). (3.11) 

n=l 

-1 - T j d\h{\) log(l + e-^) -T J dA/ 2 (A) log(l + e^) (3.12) 

;i - 2/i) - e A (A = 0). (3.13) 



The density as a function of the chemical potential is determined from ( |3.3|) . The expressions 
and ( 3. 11|) determining the free energy are equivalent to the corresponding expressions 



in as indicated in the appendix. 

3.2 Finite Temperature Excitations 

We now construct finite temperature excitations above the equilibrium state, following the tech- 
niques developed in [K|. The idea is as follows: The equilibrium state consists of an infinite 
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number of fermi seas. The state is highly degenerate, so arbitrarily choose one representative. 
Now add a particle to one of the seas, all the other seas will change in response. The energy of 
such an excitation is independent of which equilibrium state one started from, and is called a 
finite temperature excitation. 

We begin with a single particle excitation in the sea of strings of length n' . Consider an 
equilibrium state for large but finite L. If one inserts an additional particle into the sea at X p , 
this leads to an order one excitation above the equilibrium state. This is shown in Figure 1. The 
rapidities in all the seas will shift, A™ — > A™, Ap — > A.p, with A™ - A™ ~ A p - Ap ~ 0(1/ L). We 
will describe this flow by introducing shift functions [ |17] : 



A^-A 

i'n(A|A p J = lim - 

L— >oo , 

and 



S n (X\X p ) = lim (3.14) 



5 A (A|A P ) = lim , A/3 A " , (3.15) 

L^oo Ap + \ - A/3 

where the difference in the denominator is between neighboring vacancies. The energy of the 
excitation is 

oo N n 

A£ = + ££(<£&) (3-16) 

n=l o=l 

oo N n 

* 4'(Ap)-EEM(AS)(AS-A^). (3.17) 

n=l a=l 

In the thermodynamic limit the sums become integrals 

oo „ 

AE = e° n ,(X p ) -J2J dX(d x e° n )(X)9 n (X)S n (X\X p ), (3.18) 



n=l 



where 9 n (X) = (1 + e^") -1 and #a = (l + e~) _1 . In a similar way, the momentum of the excitation 
is determined to be 

A 00 r \ 

AP = P {CZ) - J2 / p'(-)e n (X)S n (X\X p ) = k n (X p ), (3.19) 

Til ~y J n 

where p'(—) = d\p(-). The shift functions can be found from the BAE. Consider the equations 
for vacancies, 

\n oo N m JV A yn A 

= 27T/" + £ £ e nm (A^ - A™) + E 0(- -) (3-20) 

n m =l 7 =l /3=1 n 
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and 

CO Wn A \n 

£ £ e(^-^SL) = 2nJ, (3.21) 



n=l ct=l 



where J" - and J now take on all integers values within the allowed range. The BAE after adding 
one particle to the n' sea are 



U(-f) = 2tt(J" + s n ) + £ e nm(AS - A™) + £ 0( g - g ) + 0„„'(A~™ - Ap) - nn ,(A2 - A,) 

(3.22) 



n m =l 7 =l p-\ n 



and 



£ *(^-^) + - 0(^^) = 2vr( J + S A ). (3.23) 



j n n n 



Here s n (n 7^ n'), s A are arbitrary half integers, s n ' is an integer arising from from a freedom to 
shift the excited state distribution by a constant. Shifting the seas by by half-integers (integer) 
leads to an 0(1/ L) shift in the energy and an 0(1) shift in the momentum. Subtracting the 
shifted equations from the unshifted one and writing the differences f(x + dx) — f(x) as f'(x)dx, 
one obtains 

00 N m 

LUK){K~K) = ££AUA£-A-)((A2-AS)-(A™-A-)) (3.24) 

m=l 7=1 

+ Efn(K-^)((K-K)-(^-~^)) (3-25) 

0=1 

- ^-®nn>CK-\p)-S n (3.26) 



and 

00 N n 

V V fJka - \ n M(\a - \») - (X" - X")) - 
n=l«=l 27F 71 

Then using ( |2.23 ), ( 2.24 ) one obtains the following integral equations for the shift functions: 



£ E /»(A/j - A^)((A, - A3) - (A™ - A™)) - _6I(-A_JV) = _ s a (3.27) 



00 . 

S n (X) = ~E duA nm (X - u)6 m (u)S m (u) (3.28) 

m=l 



dA/ n (A - A)#a(A)Sa(A) 



— nn /(A - A p ) - s r 
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and 

1 _,A-A, 



and also 



5 A (A) = I ^/„(A - \)9 n {X)S n {\) + —9{ pL) + s A (3.29) 

n=l J Zty n 

We can now show an interesting identity. First rewrite 

1 f x p 1 1 

-— 0nn'(A-A p ) = -/ di/— (0„0„„/)(A - v) - — 9 nri '(oo) (3.30) 

Z7T J-oo Z7T Z7T 

/^p 1 
rfz/A n „/(A - i/) - ~(2rnin(n,n) - 6 nn > - 1), (3.31) 
-oo 2 

f^(^-A) = - *«v(A - u) + ~ (3.32) 

27T n J-oo 2 

Taking the derivative of the integral equations for the functions e n ,e A ( (|3.7|) and (|3.8p ), with 
respect to A, A and noting that all of the kernels are symmetric functions, one can see by simple 
substitution that 

AE = e n , (A p ) - B(2n' - 1) + p. (3.33) 

In other words, the functions e n may be interpreted as the dressed excitation energy above the 
equilibrium. Proceeding in a similar way for particle-hole excitation in the A sea, one obtains 

AE = e A (A p ) - B - li (3.34) 

and 

AP = -J2 / d\p'(-)6 n (\)S n (\) = k A (A p ), (3.35) 
where the shift functions £„(A|A P ) and £ A (A|A P ) satisfy 



oo - 

S„{\) = - J2 I dvA nm (\ - v)0 m (v)S m (v) (3.36) 

m=l 



dAf n (\-A)9 A (A)S A (A) 



1 0( \- A p , 



27r n 

and 

oo „ 

S A (A) = Y, J dXf n (A - \)6 n (\)S n (\) + s A . (3.37) 



71=1 ' 
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Hole excitations may also be constructed. The equations for the shift functions are the same 
as (|3.28|) , (|3.29j) , except the inhomogeneous terms have opposite sign.. This leads to the energy 



AE = -e n ,(\ h ) + B(2n - 1) + /x. 



(3.38) 



Composite excitations are sums of the elementary ones, for instance a particle-hole excitation in 
the n' sea has shift functions 



S n (X\X p ,X h ) = S n (X\X p ) + S n (X\X h ), 
A general order one excitation above the equilibrium is of the form 

oo oo 

A£ = EE e n(A p J - EE f »W + E^a(A p J - Ee A (AO 

n=l k n=l I m n 

with momentum 

oo oo 

AP = EE U\ k ) -EE UK) + E UKJ - E ^(a^; 



(3.39) 



(3.40) 



(3.41) 



n=l k 



n=l I 



An interpretation of these results is as follows [IH]: Consider the functions e n (X,fi,T) and 
6a(A, //, T) as temperature dependent excitation energies of a noninteracting fermi system. If 
one calculates the free energy of such a system, the result is exactly the answer for the true 
free energy ( |3.11| ) of the interacting system. There is a freedom to choose how many species of 
fermions represent the system, either an infinite number for equation (|3.11| ) or two for equation 



4 Ground State and Excitations 
4.1 Ground State 

We shall obtain the ground state and excitations by taking the T — > limit of the thermodynamic 
equations. Taking the zero temperature limit of equations ( |3.7|) and (|3.8|) , we note that e n > 
for n > 2. This makes the dependence of the integral equations on e n , for n > 2, vanish in the 
zero temperature limit yielding the dressed energies 

Ci(A) = e?(A) + J dAfx(X - A)e A (A), (4.1) 
e A (A) =-{» + B) - J dXh{k - A)er(A), (4.2) 
12 



and 

e n (A) = e°(A) + J dAf n (\ - A)e A (A) - J dvf n -i{\ - i/)e7», (4.3) 

where 

.-(A)={f) ;[:<:. (4,) 

This tells us that the ground state consists of two fermi seas: one of N\ parameters and one of 
N\ parameters. This contrasts with Lai's ansatz which yields a ground state consisting of spin 
up electrons and and bound state singlet electron pairs. 

In spectral parameter space, the seas are filled according to Figure 2. The Ni sea is filled 
symmetrically around A = up to some fermi level ±Q±. The sea is is filled placing particles 
at the maximum spectral parameter down to a fermi level Q\. 

The ground state energy is 

^ = 1 - f Ql dA47rA(A)p?(A), (4.5) 



L 

where the densities pf and a p satisfy the integral equations 

rf(A) = A (A) - ( [' QA + fVA(A - A)^(A) (4.6) 

and 

<t*(A) = [ Ql dAA(A-A)p?(A). (4.7) 
The fermi boundaries Qi and Qa may be determined from 

D = -p = 1 — / dAp?(A) + ( / + / )rfAa p (A) (4.8) 

and 



dAp?(A) - ( / + / )dAa p (A). (4.9) 

Qi J-oo JQa 

In zero magnetic field, the magnetization density of the ground state is zero, due to a theorem 
by Lieb and Mattis [IB]. This requires that N\ = L — N\ — 1. The number of available spaces 
for N\ parameters according to (|2.19|) is iVi, so the sea is entirely filled. Thus the fermi sea 
boundary Qi goes to oo in the thermodynamic limit. 

Considering the half-filled case, N/L = 1, in zero magnetic field, we have Qa = and 
Qi = oo. Here the equations for the ground state are easily solved by fourier transform. The 
result for the ground state energy is E = 1 — 2 log 2, in agreement with || and ||. 
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4.2 Excitation Spectrum in Zero Magnetic Field 

We shall consider order one excitations spectrum in zero magnetic field. The dressed energy of 



the excitations is given by the functions ex, 6a satisfying (|4.1| ) and (|4.2|) . Consider the case of L 
even, Ni odd, and Aa even. The momenta of an excitation in the Ni sea is given by the zero 
temperature limit of (gT9|),(ggg),and ( ggg) : 



fcx(A p ) = p(A p ) - / )Si(HA P ) (4.10) 
where S±,Sa satisfy the equations 

+ / )dA/ 1 (A-A) ( S , A (A) ) (4.11) 

-co jQa 



and 



S A (A|A P ) = [°° dXh{A - X)S 1 {X) - -U(A - X p ). (4.12) 

Excitations in the ATa sea have momentum given by the zero temperature limit of fl3.35|) , (|3.36|) , 
and ( |Q7| ): 

fc A (A p ) = - [ Ql dXp'(X)S 1 (X\A p ), (4.13) 



where the shift functions Si, S\ satisfy 



+ / dAfi)(X - A)5 A (A|A P ) + — 9(X - A p ), (4.14) 

-oo JQ A Z7T 

S A (A\A P ) = / dA/i(A - X)Si(X\A p ). (4.15) 



For the case under considiration, all the shifts s 1 ,s A , s 1 , and s A are zero. The momentum of 
adding arbitrary numbers of particles and holes to the fermi seas is 

AP = Yl M\>)~ E h(X h )+ J2 k A (A p )- ]T fc A(A,). (4.16) 

N± particles N\ holes N\ particles N\ holes 

To see how these functions combine into excitations, first consider the less than half filled 
case, D < 1. The possible elementary processes are: 

1. Create hole in A^x sea: AN± = —1, AAa = 0. This corresponds to adding a spin up 
electron, AN e = 1 AS = ~ . The energy and momentum are 

AE = -ei(X h ), AP = -h(X h ). (4.17) 
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2. Create hole in N A sea: AN\ = 0, AN A = —1- This corresponds to removing a spin down 
electron, AN e = —1, AS" = |. This creates a hole in the jVj. sea, so the energy and 
momentum are 

AE = - £l (X h ) - e A (A h ), AP = -ki(X h ) - fc A (A fc ). (4.18) 



3. Transfer particle from N A sea to Ni sea: ANi = +1, AAa = — 1. Removing a A parameter 
increases holes in the Ai sea by one, in which the one A parameter is placed. Thus, there 
is no parametric dependence on A. This excitation removes 2 electrons, with no change in 
spin, AN e = —2, AS = 0. The energy and momenta are 

AE = -e A (A h ), AP = -k A (A h ). (4.19) 



4. Transfer particle from A"! sea to N A sea: A A"! = — 1, AN A = +1. As in type 3, there is no 
parametric A dependence. This excitation adds 2 electrons, with no change in spin. The 
energy and momentum of the excitation are 

AE = e A (A p ), AP = fc A (Ap). (4.20) 



5. Bind particles into bound state: AiVj = -2, AN A = -1, AA" 2 = +1. 

This opens up two holes in the h+ j sea, and one in the j sea. There is only one vacancy 
for the bound state, so there is no parametric dependence on its location. This excitation 
does not change the number of electrons or the spin. The energy and momentum are 

AE = -e 1 {\ hl )-e 1 {\ h2 )-e A {A h ), AP = -h(\ hl ) - h(\ h2 ) - k A (A h ). (4.21) 



4.3 Excitations at Half Filling 

In the half filled case, D = 1, we can solve the integral equations ( |4 . 1| ) , ( P~2"D , ( [4 . 1 0| ) - ( \i . 1 5| ) for the 
energy and momenta. The dressed energy and momentum of an excitation in the A^ sea are 



ex(X) 



-7T 



and 



cosh(f ) 

,7rA, 



fc x (A) = -tan~ 1 (sinh(^)) 



7T 

2' 



(4.22) 



(4.23) 
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For an excitation in the N A sea, the dressed energy and momentum are 

e A (A) =J R(A)-21og2 (4.24) 

and 

fc A (A) = / dA'R(A'), (4.25) 

J — oo 

where 

/oo p —ikw 

The n-string excitations have zero energy, e n> \ = 0. The ground state has both fermi seas filled, 
Nt = L/2, N A = L/2-1. 

These functions combine to form excitations in the following ways: 

1. Create hole in N A sea: ANi = 0, AAA = — 1. This corresponds to removing a spin down 
electron, AA" e = —1, AS = |. This creates a hole in the h+ j sea, so the energy and 
momentum are 

AE = - 6l (A A ) - e A (A h ) AP = -h(X h ) - k A (A h ). (4.27) 

2. Transfer particle from N A sea to Ni sea: = +1, AAA = — 1 In contrast with the case 
away from half filling, now two holes in the N A sea are created. This excitation removes 2 
electrons, with no change in spin, AA" e = —2, AS = 0. The energy and momenta are 

AE = -e A (A hl ) - e A (A h2 ), AP = -k A (A hl ) - k A (A h2 ). (4.28) 

3. Create holes in N\ and A^a seas: AA^ = — 1, AAA = — 1- This is not simply a sum of type 
1 and 2 excitations for D < 1 because here there is no hole in the J, sea. This excitation 
has AAg = 0,AS = 1. The energy and momentum are 

AE = -ei(A hl ) - e 1 (X h2 ), AP = -h(\ hl ) - h(X h2 ). (4.29) 

4. Bind particles into bound state: AAi = —2, AAA = — 1, AA^2 = +1. This opens up 
two holes in the Ai sea, and none in the AA sea, in contrast with type 5 excitation for 
D < 1. There is only one space for the bound state, so there is no parametric dependence 
its location. This excitation has AA" e = 0, AS = 0. The energy and momentum are 

AE = -e x (A fel ) - e^AfcJ, AP = -*i(A fcl ) - h(X h2 ). (4.30) 



The excitations 3 and 4 are identical to those in the XXX spin chain [19]. The results here 
agree with those obtained in || using Lai's and Sutherland's ansatz. 
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A Relationship to Lai's Ansatz 

In this section we shall how the solutions of Essler and Korepin's ansatz, ( |2.6|) and (|2.7| ), are 
related to solutions of Lai's ansatz (|A. 1|) and ( |A.2| ). A mathematical equivalence between the 
ansatzes has been shown in || and |20|j . 

First we recall the Lai solution 0, |5|] which treats the electrons as dynamical objects in a 
background of empty sites. The N± + iVj spectral parameters A; describe the kinetic degrees 
of freedom, the N± parameters describe the motion of the spins of the electrons. The spectal 
parameters must satisfy the BAE: 

fi(££3) / — 1, iV| + Nf, (A,, 



In this ansatz, there are two types of string solutions: 

1. bound singlet electron pairs : composed of 2 A and 1 A parameters: 

\ = A + t \' = A-i (A.3) 

2. spin bound states: take n A parameters and set 

As = A + i(n + 1 -25) 5 = l,...,n. (A.4) 

In the thermodynamic limit one define densities p,er', oi, and a n (n = 2, oo) for free electrons, 
bound pairs, down spins, and spin bound states respectively. The thermodynamic limit of the 
Lai BAE is 



//(A) = / X (A) - / dA/i(A - A)^(A) - ^ / duf n (\ - v)o*{u), (A.5) 

J n=l J 

a\A) = / 2 (A) - / dA'/ 2 (A - A')^(A') - / dA/i(A - A)p p (A), (A.6) 
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and 

ai{\) = J dvf n (X - V )fP(y) -jrj duB nm (X - !/)<(!/), (A.7) 



m=l 



where 

D _ / f\n-m\ + 2/| n _ TO |+2 + ••• + 2f n+m _2 + f n +m H n ^ 171 , 

nm ~ 1 2/ 2 + ... + 2/ 2 „_ 2 + f 2n if n = m ■ [A -* } 

Now if one makes identification between the densities in the two ansatzes 



we may see that equations ( |A.5| ), ( |A.6| ), and ( |A.7[ ) it coincide with ( |2.2b| ) and ( |2.24j ). To 
show this one only needs the identity 

fnfm — fn+m, (A. 10) 

where have written the action of the kernel as an integral operator on some function g(A): 

f n g(X)= f n {X-u)g{v). (A.ll) 



Thus the role of particles and hole solutions are interchanged in the | + | sea in Lai's ansatz 
and the N\ sea in Essler and Korepin's ansatz. One expects that this will hold in a finite box as 
well, the hole solutions of ( J2.7[ ) will coincide with the particle solutions for bound pairs in Lai's 
ansatz. 

Using this particle-hole correspondence, one can show that the free energy obtained here is 
equal to that in ref. p|. Identifying the energy functions, one can rewrite the integral equa- 



tions for the e functions, (|3.7| ) and (|3.8|) , to coincide with the corresponding equations in Lai's 
ansatz, equation (5.16) in [[J. 
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